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Abstract. - In this paper, it is shown that a configuration modulated system described by the 
Frenkel-Kontorova model can be locked at an incommensurate phase when the quantum zero 
point energy is taken into account. It is also found that the specific heat for an incommen- 
surate phase shows different parameter-dependence in sliding phase and pinning phase. These 
findings provide a possible way for experimentalists to verify the phase transition by breaking 
of analyticity. 



The Frenkel-Kontorova (FK) model was proposed to study successive commensurate- 
incommensurate phase transitions observed in many configuration modulated systems such as 
low-dimensional conductors having charge density waves and some ferroelectric materials (see 
recent review Q). It has also been used widely to study transmission in Josphson junction and 
atomic-scale friction -nanoscale tribology ^, |^ . More recently, this model has been employed to 
study transport properties of vortices in easy flow channels Q and strain- mediated interaction 
of vacancy lines in a pseudomorphic adsorbate system 

The FK model consists of a chain of particles having nearest-neighbor harmonic interaction 
and adsorbed in an external periodic potential. The Hamiltonian of the FK model is 
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which can be cast into a dimensionless one 

H = + ^(^'+1 Z^)' - 7^ cos2nX,], (2) 

where Pi and arc the dimensionless classical momentum and position of the ith particle, 
k — is a dimensionless coupling constant and /i the dimensionless lattice constant at 
vanishing k. The most important quantity of interest is the ratio (p) of the mean distance 
between neighboring particles and the period of the external periodic potential (one) 

P = J™ ^° ■ (3) 

Because the competition between two length scales the FK model exhibits abundant interesting 
phenomenal^, ^ ^ . Among many others is the mode-locking which manifests itself in the 
'devil's staircase' || I , i.e., the dependence of p on /i is described by a highly pathological 
function. Namely, p{^) has the property that for all rational number z/, there exist real 
numbers and ^{v+) such that p{pL) = v \i p,{v-) < /i < /i(z/+). These mode-locking 

intervals show up as horizontal plateaus in p{p). The widths of the plateaus are determined 
by the Farey sequence, i.e. the widest plateau between two plateaus at p — p/q and p = p' /q' 
is the plateau a,t p = {p + p')/{q + q'), where p/q and p' /q' are rational fractions. Accordingly, 
one can construct the Farey tree of the rationales with successive order starting from 0/1 and 
1/1. Therefore, assuming fc is a temperature- and pressure-dependent parameter, the model 
will exhibit a lock-in commensurate phase at relatively high order by changing temperature 
and/or pressure. This mechanism has been used to explain experimental results observed in 
thiourea and in epitaxial thin film . 

However, to obtain a reasonable and more accurate result one should consider quantum ef- 
fect, in particular in low temperature regime. This problem becomes more and more important 
when people starting working on nano systems at very low temperature such as nano-tribology. 
One may ask: What happens to the commensurate and incommensurate mode-locking when 
the quantum effect is taken into account. In other words, can devil's staircase survive the 
quantum fluctuations? This is not clear up to now. In fact, we have only limited knowledge 
coming from the continuum model, i.e. the sine-Gordon model, which is valid only for very 
small k. In the continuum limit, Bak and FukuyamaQ discussed the local stabilities of the 
commensurate phases. They found that quantum effect would destroy the normal mode-locking 
staircase if the quantum fluctuation is large enough. However, for a more general case of k the 
picture is still incomplete. 



On the other hand, the study of Borgonovi et al\ll\ by using the quantum Monte Carlo 
method shown that the quantum fluctuations, similar to the thermal fluctuations [p^ , smears 
out the discontinuity of hull function in classical FK model. One of the significant results is 
that the quantum effects, mainly from tunneling, renormalize the standard map to an effective 
sawtooth map. This phenomenon was also recovered by Herman et al [ p^ by using the method 
of mean-field theory and Hu, Li and Zhang p3[ by using the squeezed state approach. Moreover, 
Hu and Lip^ found that although the quantum effects smear out the breaking of analyticity 
transition, the remnant of this transition is still discernible in the quantum FK model, which 
is demonstrated by the crossover of the ground state wave function from an extended one to 
a localized one as the coupling constant is increased. The transition also shows up in other 
relevant parameters [ p^ . 

In this paper we concentrate on effect of quantum zero point energy. The other effects such 
as quantum tunneling will be neglected. This approximation is valid provided that Ti is not too 
large |l^. As we will see soon that the quantum zero point energy can lead to interesting 
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results. For small k, all the commensurate phases seem to be destroyed. This agrees with Bak 
and Fukuyama's finding|^]. For large k the situation becomes more complex. The plateaus 
can be destroyed and enhanced depending on the system's parameters. More interestingly, 
for suitable scale of the quantum effect |(^ the system may be locked into incommensurate 
phases. This mode-locking of incommensurate phase has a very special meaning in laboratory 
experiment. 

In the low temperature regime, we assume that the quantum effect just causes the particle 
a small fluctuation around its equilibrium position. The phonon spectrum is calculated by 
linearizing the system around its equilibrium configuration, it is determined by 

B-tj^i^O, (4) 

where B is a constant matrix with elements 

_ d^H 

and I is a unit matrix. The ground state energy is considered as the classical ground state 
energy plus the quantum zero point energy: 

f^vo + fo (5) 

where 



V - (Xf,, - Xf - uf - — ^ cos 27rXf 

,1=1 ^ ' 
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is the classical ground state energy and 
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the quantum zero point energy at the ground state. Here Xf is the equilibrium position of 
i'th particle in classical ground state. In Eq. (5) Vq is measured in unit of 6^7 while /q in unit 
of hujQ, where Wq = 7/m and h is the Plank constant. Thus the ratio h — is a measure of 
the quantum zero point energy. 

Throughout this paper, except special indication, 600 particles with periodic boundary 
condition were used in our numerical calculation. Since the system is symmetric about fj, = 0.5, 
it is sufficient to study the interval p G [0,0.5]. To find a stable phase of the system at a 
fixed /X, i.e. the phase has a minimum /, we compared three hundred phases of p = z/600, 
i = 1,2, •••,300. 

The classical ground state energy vq in Eq. (5), is determined by parameters k, p and ^i. 
It has minimums at rational values of p. p versus p shows the devil's staircase. The quantum 
zero point energy /o, however, depends only on the parameters k and p. The dependence of /o 
on p is shown in Fig. 1 for different values of fc, where one can see that /o takes maximums at 
rational values of p. It has smaller value for incommensurate phases (irrational values of p). 
Furthermore, the positions of the peaks in the figure follow the Farey sequence, i.e., between 
two peaks ai p = p/q and p = p' / q' there exists a peak at mediant p — {p + p')/{q + q') and it 
satisfies the condition min{/o(p/g),/o(p'/(?')} < M{p+p')/iq + q')] < max{/o(p/g), /o(p'/g')}^ 



(^) In the following, the "quantum effect" simply means the quantum zero point energy effect. 
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The peaks at high order rationales are indeed flattened. The orders of rationales above which 
the peaks are flattened are determined by k. And it is found that more and more peaks at 
rationales of higher orders are flattened as k is increased. The higher the order of the rational, 
the earlier the peak at the corresponding position is flattened. 

The stable phase of the system is a consequence of the competition between the two terms 
in Eq. (5), vq favorites a commensurate phase but /o an incommensurate phase. The final 
phase of the system depends on the parameters h, k, and /i. In the case of ?i ^ 1, i.e., in 
the classical limit, vq gives the major contribution, and it determines the minimum energy 
phases of the system. Thus the devil's staircase mode-locking in the classical model survive 
the quantum effect. On the other hand, in the case of relatively large h, the quantum effect 
dramatically changes the classical mode-locking structure of the system. 

In Fig. 2, we plot p, at which / has minimum, as a function of /i for different values of h at 
k — 1.1. The thick line in the figure is the result of ?7, = 0. The plot shows that those plateaus 
at ?l = evolve in two different ways as Ti is increased. Some plateaus at the rationales, whose 
corresponding peaks in the fo — p graph is flattened, may persist and even be enhanced in a 
certain interval as h is increased. For example, in Fig. 1 no peaks show up in the interval 
< p < 1/6 in curve oi k = 1.1. However, as fi increases, the plateaus at 1/7, 1/8, 1/9 and 
1/10 become more and more evident as is demonstrated in Fig. 2. Similarly, the plateaus at 
p = 2/11 and p = 3/17 expand tremendously as fi is increased, while Fig. 1 does not show any 
peaks at the corresponding positions. The expanding rate of the survived plateaus is different. 
The plateaus at rationales of higher orders expand faster. On the other hand, the plateaus, 
at the positions where fo — p graph shows peaks, are flattened quickly with the increase of 
h. These plateaus locate at the positions with rationales of lower order. The most interesting 
result is those new plateaus created at the positions of irrationals. They grow up quickly, and 
become wide plateaus in the p — p. graph, see e.g., the plateaus designated by the continued 
fraction representation of the numbers in Fig. 2. These plateaus do not exist at ?i = (thick 
line). 

According to the number theory, a number p/q can be written as a continued fraction 

- = i = [aoaia2...\. 

0,2 - ■■ 

A rational is expressed as a finite continued fraction and an irrational a nontcrminating 
continued fraction. For numerical calculation we use the Fibonacci sequence to approach 
the irrationals, for example we use p = 233/610 to approximate the golden mean [21111...] and 
p = 144/521 the silver mean [31111. ..]|l^. In these cases our computation were performed with 
610 and 521 particles, respectively. We call the plateau at an irrational number as ^an irrational 
plateau' and refer to this kind of mode-locking as incommensurate phase mode-locking in the 
following. Surprisingly, the irrational plateau is obviously related to the 'irrationality' of 
the number (in the sense that the golden mean is the most irrational 'irrational' and then 
the silver mean and copper mean etc), i.e., the more 'irrational' the number, the easier the 
system can be locked at. The physical meaning of this statement are two folds. First, by 
fixing k and increasing Ti from zero, the first locked irrational plateau appears at the golden 
mean p — [21111...] followed by the silver mean p — [31111...], and then the copper mean 
p = [4111...] and so on. The widths of the plateaus are ordered respectively, i.e., for a fixed 
fi, one has /([21111...]) > ^([31111...]) > /([41111... ])..., where l{p) denotes the width of the 
plateau at p. Second, for a fixed value of fc, the number of the irrational plateaus is finite and 
it increases as k is increased. In the case of fc > 1.5 most of the mode- locking plateaus seem 
to be located at the positions of rationals of high order. The only irrational plateau can be 
recognized is that one locating at the position of the golden mean. With the decrease of k the 
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irrational plateaus at the silver mean, the copper mean, the iron mean, . . . , and the irrational 
numbers close to these well-known irrational numbers such as [22111...], [21211...], [32111...] 
and [31211...] etc. appear gradually. At fc = 1.1 we find that the plateaus at the silver mean 
and the copper mean already shown up, as is shown in Fig. 2. Our calculations indicate 
that at fc = 0.9 an irrational plateau also appear at the position of the iron mean [51111...]. 
The irrational plateaus at [61111...], [71111...], and [81111...] show up too a,t k = 0.7, but the 
plateaus at rationals can not be seen any more when k below this value at relatively large 
h. This result, i.e. for small k all the rational plateaus are flattened as the quantum effect is 
incrased, agrees with the discovery by Bak and Fukuyama [^ qualitatively. 

The mode-locking of the incommensurate phase shown in Fig. 2 has a special meaning to 
experimentalists. It might allow us to detect the signature of the phase transition by breaking 
of analyticity||, ^ in laboratory. The idea is based on the measurement of the specific heat 
of an incommensurate phase. The specific heat of the system is given by 



where T is the system temperature measured in unit of fuvo/kB, and ks is the Boltzmann 
constant. Cy is a function of the parameter p, k and T. It does not depend on In a 
commensurate phase e.g. p = 1/3, we find that Cy decreases exponentially with k. In an 
incommensurate phase the specific heat shows a quite different behavior. It also decreases 
exponentially with fc for fc > kc{= 0.9716...). But it changes insensitively to fc for A: < kc. This 
can be easily seen from Fig. 3. Therefore, experimentally when the system is found to be in a 
locked incommensurate phase, one shall be able to observe this different parameter-dependence 
behavior in specfic heat by changing external parameters such as pressure. This might be useful 
to verify the existence of kc in real experiment and detect the manifestation of phase transition 
by breaking of analyticity. 

In summary, we have studied the effect of the quantum zero point energy on commensurate 
and incommensurate phases. If it is negligible the FK model shows the usual devil's staircase 
mode-locking, and the system can be locked in commensurate phases at relatively high order. 
In the opposite case, the scenario changes dramatically and becomes more complex. If k 
is small the mode-locking plateaus shown up in the classical model are fiattened totally by 
increasing the quantum effect. This result agrees with that of a continuum limit. However, if 
k is relatively large, the mode-locking plateaus at rationals can be flattened, persistent or even 
enhanced by the quantum effect. The occurrence of the incommensurate phase mode-locking 
depends on the degree of irrational of p. The mode-locking of incommensurate phase found in 
this paper provides us a possible way to verify the phase transition by breaking of analyticity 
experimentally. This phenomenon has not yet been confirmed in laboratory, although it has 
been predicted theoretically for two decades. We hope that our findings in this paper may 
draw attentions from experimentalists. 
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Fig. 1. - The quantum zero point energy fo as a function of p. Three different curves correspond to 
k = 0.9, f .f and f .3, respectively. The peaks happen exactly at rational numbers, while the minimums 
locate at the irrationals. The number of peaks decreases as k is increased. 



Fig. 2. - p versus p.. The p given in this figure corresponds to that value at which / takes minimum. 
The curves are ordered from right to left on the left hand side of point P while from right to left on 
the right hand side of the point as h is increased from 0, 0.05, 0.1, 0.2, 0.3, 0.4. 
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Fig. 3. - Specific heat C„ versus k for an incommensurate phase with the golden mean p. The 
temperature is fixed at T = 0.2 and p — 233/610. The curve is obviously divided into two parts. 
At A; < fee = 0.9716... region, Cv is marginally independent of fe, while at fc > fee the Cv decreases 
exponentially. 
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